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Abstract. The tails of gravitational waves are caused by scattering of linear waves 
onto the space-time curvature generated by the total mass-energy of the source. Quite 
naturally, the tails of tails are caused by curvature scattering of the tails of waves 
themselves. The tails of tails are associated with the cubic non-linear interaction 
between two mass monopole moments and, dominantly, the mass quadrupole of the 
source. In this paper we determine the radiation field at large distances from the 
source for this particular monopole-monopole-quadrupole interaction. We find that 
the tails of tails appear at the third post-Newtonian (3PN) order beyond the usual 
quadrupole radiation. Motivated by the need of accurate templates to be used in the 
data analysis of future detectors of gravitational waves, we compute the contribution 
of tails, and of tails of tails, up to the 3.5PN order in the energy flux generated by 
inspiraling compact binaries. 
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1. Introduction 



1.1. Motivation and overview 

Gravitational waves propagating through vacuum from their source to infinity 
share all possible contributions associated with products of multipole moments - 
indeed, this is a consequence of the infinite non-linearity of the field equations. 

At the quadratic non-linear order, two multipole interactions play a prominent role 
(and yield effects which are representative of that order). The first interaction is that 
of the (mass- type) quadrupole moment M pq (t), which dominates the radiation field for 
slowly-moving sources, with itself. See [1] and references therein for the computation 
and discussion of this interaction (henceforth we shall refer to [1] as paper I). The 
second interaction is between M pq (t) and the static mass monopole moment M of 
the source (Schwarzschild mass, or, more precisely, ADM mass of the source). Such 
interaction M x M pq (t) is physically due to the propagation of quadrupole waves on the 
Schwarzschild background associated with M. In particular the scattering of waves 
onto the potential barrier of the Schwarzschild metric produces the so-called tails, 
which are pieces of the field depending on the parameters of the source at all instants 
from — oo in the past up to the retarded time t — r/c. 

The numerous works related or devoted to tails include some mathematical 
investigations of curved space-time wave equations [2-9], several investigations and 
constructions of post-Minkowskian expansions [10-17], the studies of the linear 
perturbations of the Schwarzschild metric by fields of various spins [18-24], some 
discussions of physical properties of tails [25-31,53], and the development of accurate 
wave-generation formalisms [54,55,32-35]. Not only the tails exist as theoretical 
objects predicted by general relativity, but they should exist in the future as real 
observed phenomena. Indeed the presence of tails in the gravitational-wave signals 
generated by inspiraling compact binaries should be deciphered by the planned 
experiments VIRGO and LIGO (see [36,37,30,31]). 

In the present paper we develop further this subject by computing the cubic 
interaction between the quadrupole M pq (t) and two monopoles M. Physically this 
"monopole-monopole-quadrupole" interaction is responsible for the scattering of the 
linear quadrupole waves M pq onto the second-order (M 2 ) potential barrier of the 
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Schwarzschild metric, and for the scattering of the quadratic tails M x M pq themselves 
onto the first-order (M) potential barrier. The latter effect produces the so-called 
"tails of tails" of gravitational waves, as they can pictorially be referred to. In fact we 
shall employ this crude appellation for the whole M 2 x M pq interaction. 

Like the quadratic tails, the cubic tails of tails could be computed using black- 
hole perturbation techniques. In this paper, we rather employ the particular post- 
Minkowskian approximation method proposed in [16], and which is used to compute 
the quadratic metric M pq x M rs in paper I. In Section 2 we recall from previous 
work some relevant results on linear and quadratic metrics, and we determine the 
cubic source term in the field equations (in vacuum) corresponding to the looked-for 
interaction M 2 x M pq . Then we compute in Section 3 the (finite part of the) retarded 
integral of this cubic source term, restricting ourselves to the leading order 1 jr in the 
distance to the source. The complete radiation field for the M 2 x M pq interaction is 
obtained, and discussed, in Section 4. 

We find that the tails of tails carry a supplementary factor 1/c 6 with respect to 
the dominant quadrupole radiation, and therefore contribute to the radiation field at 
the so-called third post- Newtonian (3PN) order. Such a high post-Newtonian order is 
a priori quite small in absolute magnitude, but it is still relevant to the observations 
of inspiraling compact binaries by VIRGO and LIGO [36-40,58]. Essentially the post- 
Newtonian corrections in the field affect through gravitational radiation reaction the 
evolution of the binary's orbital phase, and the latter observable will be monitored 
very accurately in future detectors thanks to the large number of observed periods 
of rotation. In the author's opinion, it is remarkable that such a cubically non-linear 
effect as a tail of tail should be known in advance for comparison with real observations 
- and, therefore, should in principle be detected at the same time. In Section 5 we 
compute the tails of tails at 3PN order occuring in the total energy flux generated by 
inspiraling compact binaries (the energy flux is the crucial quantity to predict because 
it yields via an energy balance argument the effects of radiation reaction) . We compute 
also the contribution of quadratic tails at the 3.5PN order (extending previous results 
at the 1.5PN and 2.5PN orders). Besides the tails of tails at 3PN in the energy flux of 
binaries, there are also some "instantaneous" contributions, the computation of which 
is in progress [56,57]. 

In the particular case where the mass of one body is very small as compared with 
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the other mass, the radiation field of compact binaries has been computed analytically 
up to a very high post-Newtonian approximation [48-51]: notably up to the 4PN 
order [50] and, more recently, 5.5PN order [51]. Our results in Section 5 are in perfect 
agreement, wherever the comparison can be made, with the latter works. 

The notation and conventions are essentially the same as in paper I (a short 
summary is provided in [41]). In order to reduce clutter we pose G = c = 1 when 
indicating the G's and c's is not essential. For a short review on the post-Minkowskian 
method we refer to Section 2 of paper I. 



1.2. Notation for the Geld equations 

Our basic field variable is h af3 = y^gg af3 — r] at3 \ with g af3 the contravariant 
metric, g the determinant of the covariant metric, and rj a ^ the Minkowski metric in 
Minkowskian coordinates (signature — h ++). Subject to the condition of harmonic 
coordinates, 

d p h ap = o , (i.i) 

the vacuum field equations read 

Uh aP = N ap (h, h) + M ap (h, h, h) + 0(/i 4 ) , (1.2) 

where □ denotes the flat d'Alembertian operator, and where the source term in the 
right side is made of an infinite sum of quadratic, cubic, and so on, functionals of h al3 
and its first and second derivatives. The quadratic and cubic terms are given by 



N ap {h, h)=- h^d^h af3 + \d a h^h" v - -d a hd p h 



4 o 2, 



(1.3) 



M ali {h, h, h) = -h^(d a h^ p d l3 K + d P h«d p hP - dutfdvhP") 

+ h af3 -\d p h pv d p W v + \d p hd»h + \d p h vp d v W p 

4 o A 
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+ v 



a/3 



(1.4) 



The higher-order terms are even more complicated but will not be needed in this 
paper. 



2. The monopole-monopole-quadrupole source term 



2.1. The linear and quadratic metrics 

We look for a solution of the equations (1.1)-(1.4) in the form of a post-linear (or 
post-Minkowskian) expansion 

h af3 = Gh af3 + G 2 h af3 + + ^ ^ ^ 

where G denotes the Newton constant. The harmonic coordinate condition (1.1) 
implies that all the coefficients of the G n, s are divergenceless. On the other hand, 
the field equations (1.2) imply that the coefficients of any G n, s obey a d'Alembertian 
equation whose source is known from the previous coefficients, i.e. coefficients of the 
G m, s where 1 < m < n — 1 (see e.g. Section 2 in paper I). 

Our starting point is the linearized metric defined by the equation (2.3) of 
paper I. This metric is in the form of a multipolar series parametrized by symmetric 
and trace-free (STF) mass-type multipole moments Ml {& > 0) and current-type ones 
Sl > 1) [41]. These moments reduce, in the Newtonian limit c — > oo, to the 
usual Newtonian multipole moments [14] . As we are ultimately interested only in the 
cubic interaction between two monopoles M and the quadrupole M pq , we retain in 
the linearized metric only the terms involving M and M pq . Accordingly we denote 
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The monopole term reads 



, a/3 _ , a/3 , , a/3 2 n 



>*(M) = -4r- x M , (2.3a) 

l (M) — (M) 



^ = ^(m\ = • (2-36) 



This is simply the linearized piece of the Schwarzschild metric in harmonic coordinates, 
for which only the 00 component of our field variable is non-zero. The quadrupole 
term in (2.2) is 



h° { ° Mpq) = -2d ab [r-'Matit - r)} , 



l,0i _ 90 
n (M pq ) - L °a 



(2.4a) 
(2.46) 
(2.4c) 



We use the same notation as in (2.3) of paper I, which gives the complete linearized 
metric including all multipole terms. In the following we need rather the metric (2.4) 
in expanded form, where the spatial derivatives acting on both r _1 and t — r are 



worked out. We have 



h° { M pq ) = -2n ab r-* {3M a6 + 3rM« + r 2 M%> } , (2.5a) 
hf Mpq ) = -2n a r" 2 {m£> + rM<J>} , (2.56) 
h% pq) = -2r-i M V . (2.5c) 

Henceforth we generally do not indicate the dependence of the moments on t — r. 

Consider the quadratic metric generated by the linear metric (2.2)-(2.5). It is 
clear that involves a term proportional to M 2 , the mixed term corresponding to 
the interaction M x M pq , and the term corresponding to the self-interaction of M pq . 
Thus, 

, a/3 _ , a/3 . , a/3 , , a/3 /r, fi \ 

n 2 ~ + a {MM pq ) + a (M pq M rs ) ■ 



The first term is the quadratic piece of the Schwarzschild metric in harmonic 
coordinates, 

fc(£f») = "7r" 2 M 2 , (2.7a) 

hf M2) = , (2.76) 

h\ J M2) = -n ijr - 2 M 2 , (2.7c) 

with four-dimensional trace (h = riaph^ 13 ) 

h {M2) = 6r" 2 M 2 . (2.7d) 

The term y which constitutes the dominant non-static multipole interaction 

at the quadratic order, obeys a d'Alembertian equation whose source is given by (1.3) 
where (2.3) and (2.5) are inserted, and which can be written, with obvious notation, as 
N a P(h( M ), h( Mpq )) + N af3 (h( Mpq ), h(M))- The solution of this d'Alembertian equation, 
and, then, the complete metric ^ itself, is obtained using the present method. 

The monopole-quadrupole metric is the analogue of the 2-2 metric in Bonnor's double 
series method [10]. We simply report the result of its computation, which can be 
found in Appendix B of [32]: 

M" 1 ^) = n ab r~ A {-21M a6 - 21rM# + 7r 2 M$ + 10r 3 M^} 

/+oo 
dxQ 2 (x)M^>(t-rx) , (2.8a) 

M-'hf MMpq) = n iab r-* {-M<? - rMif - \r 2 M^ 

+ n a r- 3 {-5m£> - brM^f + yr 2 M<? j 

/ + oo 
dxQ^M^it-rx) , (2.86) 

M-^ MMpq) = n ljab r~* {-^M a& - ^rM« - 3r 2 M<? - \r*M% 

+ W" 4 [-\ Mah - \rM$ - 2r 2 M$ - 
+ n a(i r" 4 [6M j)a + 6rM« + 6r 2 M$ + 4r 3 M$} 
+ r" 4 | -M^ - rM^f - 4r 2 M^ 2) - yr 3 M^ 3) 



/+oo 
dxQ (x)M$\t-rx) , (2.8c) 

with four-dimensional trace 

M~ x h {MMvq) = n ab r~ 4 [l8M ab + lSrM^ - 10r 2 M%> - 12r 3 M®>} 

/+oo 
dxQ 2 {x)M<£{t-rx) . {2 Ad) 

The metric is composed of two types of terms, instantaneous terms depending on 
the quadrupole moment at time t — r only, and non-local (or hereditary) integrals 
depending on all instants from — oo in the past to t — r. The usual tail effects are 
contained in the non-local integrals of (2.8). Note that the non-local integrals come 
exclusively from the source terms whose radial dependence is r~ 2 [see Section 3 of 
paper I and (3.1) below]. The integrals are expressed in (2.8) by means of the Legendre 
function of the second kind Qi (with branch cut from — oo to 1), which is related to 
the Legendre polynomial Pi by 

Qi{x) = \j_ 1 Pi{v) x^y (2 - 9fl) 
= \Pt{x)\n - £ W^P^x) (2.96) 

(the first of these relations being known as Neumann's formula for the Legendre 
function, see e.g. [52]). See also (A. 15) in Appendix A for still another expression of 
the Legendre function. For future reference we quote here the expansion of Qi when 
x — > 1 (with x > 1), 

Q e (x) = -\ In ( X -^-\ - E 7 + °K X ~ !) ln (* " !)1 • ( 2 - 9c ) 

V J 3=1 3 

[On the other hand, recall that the Legendre function behaves like l/x i+1 when 
x — > +oo.] Finally the term M ^ in (2.6) is the quadrupole-quadrupole metric 
whose computation has been the subject of paper I (we do not need this term in the 
present paper). 

Now the cubic metric is made out of all possible interactions of three 
moments chosen (with repetition) among M and the quadrupole M pq , and is therefore 
constituted of four terms, 
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a 3 - a {M-i) + n (M^M pq ) + a (MM pq M rs ) + n (M pq M rs M tu ) • ^ Z - iU J 

Of these terms only the first one is known within the present approach (before this 
paper): this is the cubic piece of the Schwarzschild metric that we give here for 
completeness, 

fc(M3) = -8r" 3 M 3 , (2.11a) 
= ^3, = . (2.116) 

The term h^ 2M ) is the monopole-monopole-quadrupole metric which will be dealt 
with in the present paper. The two last terms, which involve at least the interaction 
of two quadrupole moments, will be left undetermined for the time being. 



2.2. Expression of the cubic source 

The metric h^ 2M ) obeys the harmonic-coordinates condition dph^ 2M ^ = 
and a d'Alembertian equation whose source exhausts all possibilities of generating the 
multipole interaction M 2 x M pq by means of linear and quadratic metrics. From (1.2) 
one has 

n /a/3 _ A a/3 ,cy 12 s 

LJ,L (M 2 M pq ) ~ iv (M 2 M P9 ) ' l z - lz J 

where the source term A"^ 2M ^ is obtained from N af3 and M Q/3 defined in (1.3) and 
(1.4) as 



A< (MiM pq ) = ^ (V 2 )'V P9 )) + N " P i h (M pq ), h { M*)) 

+ iV a/3 (/l( M ), /l(M„,)) + N af3 (/l(MM M ))V)) 

+ M Q/? (/i(M), fr(M), h( Mpq )) + A^ Q/3 {h(M),h( Mpq ), h( M )) 

+ M^(h (Mpq) ,h (M) ,h {M) ) . (2.13) 

In the first and second lines a linear metric is coupled to a quadratic one, while in 
the third and fourth lines three linear metrics are coupled together. The metrics have 
been distributed with all possibilities on each slots of the non-linear sources (1.3) and 
(1.4). 
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Using the explicit formulas (2.3), (2.5), (2.7) and (2.8), one obtains the cubic source 
term A^. 2M ^ after a tedious but straightforward computation. As the monopole- 
quadrupole metric (2.8) involves some non-local integrals, so does A?^ 9 , nj \ which 



can thus be split into a local (instantaneous) part, say M y and a non-local (tail) 
part, T^ M 2 M ^y. 



L (M 2 M pq ) 



a(3 



= J"' 8 



1Y (M 2 M pq ) J (M 2 M pg ) 

The result for the instantaneous part is 



a/3 



(M 2 M pq ) ■ 



(2.14) 



M -2 rOO 

M L {M 2 M pq ) 



n ao r" 7 j-516M a6 - 516rM^ - 304r 2 M^ ) 

- 76r 3 M n (3 1 ) + 108r 4 M^ } + 40r 5 M o ( ^ I , 



m ~ 2 r 0i 



a6 ^ a6 

= n iab r- 6 {4M« + 4rM<J> - 16r 2 M^ + ^r 3 M^ 



(2.15a) 



-r 4 M 



at 



, ; i 372 „ r (u 372 , . 

" ' M ( ) rM K 

5 5 



(2 > zi!z,-' w ,:!: 



84 



124 



(2.156) 



/lf~ 2 7" ij 

JW 1 (M 2 M pq ) 



92 



,3^(5) 



-» i 160 , r (2) 176 , rs) 596 9 , r^) 160 o , 



+ n a{i r 
+ r 



-5 



312 248 _ (3 ) 400 . (4 ) 104 o (5) 
M\ J rM\ j H r 2 M\ j H r 3 M\ j 



•'J -12M<?> - — rM.( 3) - -r 2 M (4) - -r 3 M.^ 



(2.15c) 



15 13 5 " ~ lJ 5 

(we recall e.g. that hij a b denotes the STF projection of ni jao = ninjn a rib [41]). 
The tail part is composed of sums of products of local terms with derivatives of 
tail terms. Using some elementary properties of the Legendre function [namely 
xQi(x) = §Q 2 (^) + lQo(x) and xQ 2 (x) = ^Qs(x) + %Qi(x)], we obtain 



M~ 2 T™ 2Mpq) = n ab r~ 3 



C dx { 



96QoM$> + 



272^ 168^ 
-^Qi + ^Qs 



rM 



(5) 
ab 



+ 32Q 2 r z M t 



■2 7i /r (6) 
ab 



(2.16a) 
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M 2 T®l I 2 Mvq } — fliabT 3 



+ n a r~ 



+00 



1 

+00 



dx{ -32Q 1 M { a t ) + 



32 8^ 



rM 



(5) 



+ 32Q 1 r 2 M^ ) 



192^ 112 _ 



rM 



ab 
(5) 



(2.166) 



M~ 2 T ij 

1V1 1 {M^M pq ) 



— flijabT I 
Jl 



+00 



dx<[ -32Q 2 M^ + 



5 5 



rM 



(5) 
ab 



/+00 r 22 

cfe<| -—Q 2 M i J ) ) + 



7 



208^ 24^ 
— <5i + yQ 3 



rM 



(5) 



/•+00 



2112 192 



rM 



(5) 



+ 



+00 



^|q 2 m^ + 



1536 o 96 o 



rM, 



(5) 
«7 



+ 32Q r 2 Mf } 



(2.16c) 



(where the Legendre functions are computed at x and the moments at t — rx). At 
this stage we have a good check that the computation is going well. Indeed A"^. 2M ^ 
is the source of the third-order field equations in harmonic coordinates and therefore 
should be divergenceless. The divergence of the tail piece (2.16) is computed using 
the Legendre equation (1 — x 2 )Q"i(x) — 2xQ' e (x) +£(£+ l)Qt{x) = S+(x — 1), where 
the distribution 5 + is defined by J^~°° dx5 + (x — l)0(x) = 0(1) for any test function 0. 
We find that all the tail integrals disappear, and get 



M ~% T (M^ Mpq) = W 4 {-48MW - 48rM<? - 16r 2 M^} , (2.17a) 
M ~ 2 9pT^ M2Mpq) = n !fli r- 4 {-8M;f - 16rM<?} - 32n a r" 2 Mif . (2.176) 

On the other hand the divergence of the instantaneous piece is computed directly from 
(2.15). We readily obtain that it cancels exactly (2.17) so that the required condition 

9 P A (Sn Mpq ) = (2-18) 

is fulfilled. 
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3. The retarded integral of the source term 

The previous check being done we confidently tackle the difficult part of the 
analysis, namely to find the inversion of the d'Alembert equation with source term 
given by (2.14)-(2.16). We shall limit ourselves to the computation of the metric at 
large distances from the source (r — > oo with t — r =const), keeping only the dominant 
1/r term at infinity (possibly multiplied by some powers of lnr). This is sufficient in 
view of applications to astrophysical sources of gravitational radiation. 

3.1. Integrating the instantaneous terms 

The general form of the terms composing (2.15) is hir~ k F(t — r), where hi is 
equivalent to a spherical harmonics of order £, and where the radial dependence is 
such that k > 2. So let us review the formulas required to compute the 1/r (and 
In r/r) terms at infinity of the retarded integral of any source term hir~ k F(t — r). 
Recall that in general the retarded integral cannot be applied directly on hir~ k F{t—r) 
because of the singular behaviour when r — > 0. We follow the procedure proposed in 
[16] to obtain a particular retarded solution (also singular when r — > 0) of the wave 
equation. It consists of multiplying first the source term by a factor (r/r ) B , where B 
is a complex number and ro a constant length scale, thereby defining a S-dependent 
fictituous source which, for large values of the real part of B, is regular, and in fact 
tends to zero, when r — > (i.e. the singularity at r = is killed). For large values 
of Real(-B) one is thus allowed to apply the retarded integral on the fictituous source 
(there is no problem at the bound r — > oo of the integral because the moments are 
assumed to be constant in the remote past). In this way one defines a function of B, a 
priori only for Real(S) large enough, but the point is that this function is extendible 
by analytic continuation to all complex values of B, except at some integer values 
(including in general the value of interest B = 0) , where it admits a Laurent expansion 
with some poles. Now it has been shown [16] that near the value B = the finite 
part of the Laurent expansion is a particular solution of the d'Alembertian equation 
we wanted to solve. We call the operator giving this solution the finite part of the 
retarded integral, and denote it by FP s = oD^ 1 (r /ro) B , or more simply by FPD" 1 . The 
finite part procedure is especially convenient when doing practical computations. See 
Appendix A in paper I for a compendium of formulas, obtained with this procedure, 
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enabling the computation of the quadratic non-linearities. More complicated formulas 
to compute the cubic non-linearities are reported in Appendix A of the present paper. 

For a source term of the type fiLr~ k F(t — r) three cases must be distinguished (as 
we have already k > 2): k = 2, 3 < k < £ + 2, and £ + 3 < k (see Appendix A of 
paper I). The case k = 2 corresponds to a retarded integral which is convergent (so 
the finite part at B — is unnecessary), and given by a non-local integral admitting 
an expansion when r — > oo, t — r = const in powers of 1/r with a logarithm of r. The 
formula, already used in (2.8) to express the tail integrals, reads 

/+oo 
dxQ e (x)F(t-rx) , (3.1) 

where Qi is the Legendre function (2.9). The leading term at infinity is obtained by 
using the expansion of Qi as given by (2.9c). We obtain [32] 



□5' [n L r^F(t - r)] = % I drF(t - r - r) [in (£) + ^ 1 



+ o 



lnr 



(3.2) 



In the case 3 < k < £ + 2 (thus I > 1), we find that the independent integral is finite 
(no pole at B = 0), and is given by a simple local expression, without logarithms. The 
1/r term at infinity is 



[(rlr»)°n L r-»F(t - r)] ,„ = - ^^S^^tt - r) 

+ 0(1) . (3.3) 

In the last case k > £ + 3, the S-dependent retarded integral admits truly a polar part, 
and its finite part is given like in (3.2) by a non-local integral, but now the expansion 
at infinity involves no logarithms of r (instead it involves the logarithm of the constant 
ro )• We have, for the 1/r term, 



FPD- 1 [n L r- k F(t - r)] 



( _ )fc+ , 2fc _3 (fc _ 3) , ^ 

(k + £-2)!(k-£-3)l r 
14 




With the formulas (3.2)- (3.4) it is straightforward to obtain the 1/r and In r/r 
terms in the finite part of the retarded integral [in short FPD^ 1 ] of (2.15): 




where the moments are evaluated at time t — r — r. Recall that in the case of the 
quadratic non-linearities, the non-local integrals come only from the source terms 
having k = 2 (see Section 3 of paper I). This is not true in the case of the cubic 
(and higher) non-linearities: besides the integrals generated by the terms k = 2, there 
are some non-local integrals generated by terms having k > I + 3. These integrals, 
given by (3.4), depend on the constant length scale tq. Thus the particular solution 
which is picked up by the finite part procedure depends on the constant ro used in 
its definition. This is quite normal, and, of course, not a problem - one needs only 
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to be consistent in using this particular solution, notably when relating the multipole 
moments which parametrize it to the source variables. 



3.2. Integrating the tail terms 

The effects which are physically associated with tails of tails come from the 
retarded integration of the tail part of the cubic source, that is T^ 2M ) given by 
(2.16). The problem is to find the (finite part of the) retarded integral of a source 
term involving a non-local integral with some Legendre function Q m , namely 



/+oo 
dxQ m (x)F(t - rx) 



(3.6) 



We study only the needed cases which correspond to k > 1 (note that with this 
notation the actual radial dependence of the source when r — > oo, t — r = const is 
l/r fc+1 ). Of course the problem is more complicated than in Section 3.1 because the 
source term at a given value of t — r is a more complicated function of r. 

The detailed computation of k,m^ l is relegated in Appendix A. Here we report the 
main results. The case k = 1 is the most interesting because it leads to qualitatively 
new results with respect to the quadratic non-linear order. When k = 1 we find the 
solution in closed analytic form, 



/+oo 
dyF^it-ry) 

x ^(y) J\xQ m (x)^(x) + P £ (y) dxQ m (x)^(x)^ , (3.7) 

where the time anti-derivative is defined by F( _1 )(t) = dt'F(t') (all the functions 
involved are zero in the remote past by our assumption of initial stationarity) . The 
solution (3.7) has been obtained thanks in particular to the mathematical formula 
(A. 5) in Appendix A. To leading order when r — > oo, t — r = const the second integral 
in the brackets dominates the first one. Inspection of (2.16) shows that when k = 1 
we need only the case m = £. Because the second integral in (3.7) can be explicitly 
worked out in this case, we have a good simplification, since notably the dominant 
term at infinity follows simply from the expansion (2.9c) of the Legendre function. 
We find a non-local expression with lnr and ln 2 r terms: 
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,^ L = -^£ 00 drF(- 1 \t-r-T) 



x 



l n 2 



+ o(r £ - 2 ) . (3.8) 



For the remainder we use the notation o(r e ~ 2 ) with a small o-symbol to mean that 
the product of this remainder with the factor r 2_e , where e is such that < e <C 1, 
tends to zero when r — > oo. This notation is simply to account for the presence of 
powers of logarithms in the expansion at infinity. 

In the case 2 < k < £ + 2, the result is simpler (see Appendix A) in the sense that 
it is given by a local expression to order 1/r, 



k ^L = - k, m a £ ^F<*" 3 >(t -r) + o (r £ " 2 ) 



However, the coefficient is still complicated: 



/+oo /-+oo f, _ x \k-3 

dxQ m (x) J dz jj—^y^ Qi{?) ■ ( 3 - 10 ) 

The numerical values of the fc, m a/s for general k, £, and m are computed in (A.17)- 
(A.18). Here we need only some values corresponding to k = 2 and k = 3 [see (2.16)]. 
These are presented in Tables 1 and 2 respectively. 

In the last case k > £ + 3, we obtain again a non-local integral, but with simply a 
1/r term without lnr (instead, there is a lnr ), 



(3.9) 



k,m^L = -^ ! °° drF^\t-r-r) 
r Jo 

The coefficients k,mPe an d k,mlt are also rather involved, 



+o(r £ - 2 ) . (3.11) 



k.rah 



k.mll 



dxQ m (x) / dz— — — Pt(z) , 



-1 



(k-3) 



i r + °° r 1 (z-x) k ~ 3 

= 2 J, dxQMjj^-j^rP^) 



In 



fc-3 1 i 
X — Z \ 1 



(3.12a) 
. (3.126) 



Fortunately these coefficients are needed only for the particular set of values k = 3, 
m = 2 and £ = 0, corresponding to the only term present in (2.16) in this category, 
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which is the antepenultimate term in (2.16c). From the formulas (A. 22) of Appendix 
A we find 



3,2 



(3.13) 



1 7 

6 ' 3 ' 27 ° ~ 72 ' 

With the expression (3.8), the values of coefficients in Tables 1 and 2, and the two 
values (3.13), we have the material for the computation of the retarded integral of 



nCtf3 

(M 2 M pq ) 



. The result reads 
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1 r23 
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(4) 



52 



r(4) 



30 



+ o (r £ " 2 ) 



(3.14c) 



All the logarithms (and logarithms square) are computed with r/2r except for one, in 
the last equation (3.14c), which is computed with r/2r and is issued from the formula 
(3.11). The other logarithms, and logarithms square, are issued from (3.8). 



4. The monopole-monopole-quadrupole metric 



4.1. The metric in the far zone 

The retarded integral of the monopole-monopole-quadrupole source, namely 
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a/3 

U {M^M pq ) 



^rKm^mJ > (4-1) 

is obtained simply as the sum of (3.5) and (3.14). Let us now follow the method 
proposed in [16] (see also Section 2 of paper I), and add to u"j^ 2M ^ a supplementary 



term v^ 2M ) so designed as (i) to be a solution of the homogeneous wave equation, 



and (ii) to be such that the sum of u"j^ 2M ) and v^ 2M ) is divergenceless. With (i) 
and (ii) satisfied, a particular solution of the cubic-order field equations in harmonic 
coordinates is u^ 2M ^ +v^fj2 M ) • Actually we do not follow exactly the construction 
proposed in [16], but adopt the slighly modified construction of 



a/3 



\M 2 M pg ) 



defined by the 



equations (2.11) and (2.12) of paper I. Furthermore, as we control only the dominant 



behaviour at infinity of the term u^ 2M ^ , 



we must check that this weaker information 



still permits the construction from the divergence of u^ 2M ) of the corresponding 



1/r term in v^ 2M y This poses no problem, and the relevant formulas can be found 
in Appendix B. The divergence of u^ 2M y computed from (3.5) and (3.14), reads 



dp 
dp 



1V1 U {M*M pq ) 
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(4.2a) 



600 J 



(4.26) 



Then the formulas (B.2)-(B.5) in Appendix B give the 1/r term in v"^ 2M ^ as 



M -2 00 

M v {M^M pq ) 

M v (M*M pq ) 



O i 1 



105 r a * V^V 



iW V (M*M pq ) 
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The complete monopole-monopole-quadrupole metric, defined by 



(4.3o) 
(4.36) 



(4.3c) 



n (M 2 M pq ) - U (M 2 M pq ) ^ U (M 2 M pq ) ' 



(4.4) 
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is therefore obtained by adding up the expressions (3.5), (3.14) and (4.3). We obtain 

^|-41n 2 (^)-41n(^-) 



M~ 2 U 



M~ 2 h 




(the moments in the integrands are evaluated at t — r — r). 

4.2. The observable quadrupole moment 

The computation we have done so far uses harmonic coordinates, which are 
convenient for constructing solutions by means of a post-Minkowskian algorithm 
(essentially because all the components of the field obey some wave equations). 
However, the harmonic coordinates entail a small disadvantage, namely the associated 
coordinate cones t — r deviate by a logarithm of r (in first approximation) from the 
true light cones along which gravity propagates. As a result, the expansion of the 
metric when r — > oo, t — r = const involves besides the normal powers of 1/r some 
powers of the logarithm of r [42,43,16,17]. This is clear from the previous result 
(4.5). The logarithms can be gauged away by going to some radiative coordinates 
which are such that the associated coordinate cones T — R (where R = |X|) become 
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asymptotically tangent to the true light cones at infinity. In this paper it will be 
sufficient to check that the coordinate transformation 

X» = + Gg* , (4.6a) 
where x M are the harmonic coordinates and £ M is given by 

£° = -2Mln(-^) , f = 0, (4.66) 

does remove all the logarithms of r in the particular case of the interaction M 2 x M pq 
(and at leading order at infinity). Note that we have introduced in the coordinate 
transformation (4.6) the same constant tq as used in the definition of the finite part 
process in Section 3. Actually we could have introduced any constant T\. However the 
choice r\ = ro, which is simply a choice of gauge (equivalent to a choice of the origin of 
time in the far zone), is especially convenient, as it will simplify some formulas below. 
Under the coordinate transformation (4.6) the metric is changed to 

where we keep only the terms corresponding to the interaction M 2 x M pq and neglect 
all sub-dominant terms o(R e ~ 2 ). Both sides of (4.7) are expressed with the radiative 
coordinates . We substitute in the right side the harmonic-coordinates linear metric 
(2.5), quadratic one (2.8), and cubic (4.5), and find that all logarithms disappear to 
order 1/-R, so that we obtain the "radiative" metric 

M-^ TM , Mn) = ^ f (Tr - r) {-4 (£) + gln (^) 

m-XW„ } = ¥ C (r « - T) { - ln fe) - m } 



71361 
2205 J 
(4.8a) 



227241 
7875 J 
(4.86) 
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+ I f~ (T B - r) ( -4 .n^ V ± in ( ^ 198 



^io ' J I V 2 ^o/ 15 V 2r oy 35 

+ o(R £ - 2 ) . (4.8c) 

The retarded time is denoted by Tr = T — R and the direction to the observer by 
N a = N a = X a /R. Note that 

K » H (m>m pq) = o (R £ - 2 ) , (4-9) 
where K v = (— 1, iV l ) is the Minkowskian null direction to the observer. 

From the radiative metric (4.8) we extract the "observable" multipole moments 
which are the quantities measured by an observer located at infinity. The observable 
multipole moments Ul and Vl parametrize the algebraic transverse-tracefree (TT) 
projection of the spatial metric in radiative coordinates, 

(H^) TT = - -V ijab Y^p[NL-2U ijL - 2 (T R ) - — Ar aL _ 2£a6(i y j>L _ 2 ( Ti j 

+ o(±), (4.10) 



R 2 ) 

where the TT projection operator is 



V ljab (N) = (Sta - NiN a )(6 jb - N.m) - htfij - NiNjWat - N a N b ) . (4.11) 

[As we use the field variable H lJ instead of the covariant metric, the formula (4.10) 
differs by a sign from (5.1) in paper I.] The moments Ul and Vl agree at the linearized 
order with the £th time-derivatives of the moments Ml and Sl [14]. 

Working out the TT projection and comparing the result with (4.10), we 
readily find that the monopole-monopole-quadrupole metric (4.8) contributes to the 
quadrupole observable moment Uij, and only to this moment, by the expression 

SU ii(TR) - ^ f «tf {TR - T) fa) + I ln (S.) + ^} . 

(4.12) 
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We add back the necessary powers of G and 1/c [recall that 5Uij has the dimension of 
(mass) (length/time) 2 ], and find that the term (4.12) carries in front a factor G 2 /c 6 , 
and therefore represents a small modification of the lowest-order quadrupole radiation 
at the level of the third post-Newtonian (3PN) order. Let us prove that there is no 
other contributions in Uij at this level. We know from dimensional arguments that a 
non-linear term in Uij involves necessarily a factor i/ c 3 ( n - 1 )+ E ^i+ s - 2 j where n > 2 is 
the order of non-linearity, where S£j is the sum of multipolarities of the n moments 
Ml and/or Sl composing the term (1 < i < n), and where s is the number of current- 
type moments Sl (see for instance Section V in [35]). Furthermore, + s — 2 = 2k 
where k is the number of contractions of indices among all the indices carried by the 
moments. For a term at the 3PN order we thus have 3(n — 1) + 2k = 6, which has 
the unique solution n = 3 (cubic interaction) and k = (no contraction of indices). 
Then the multipolarities satisfy £\ + £ 2 + £3 + s — 2 = 0, so we have necessarily 
(^1,^2, ^3) = (0,0,2) and s = 0, which corresponds indeed to the sole interaction 
M 2 x Mij. 

The equation (5.10) of paper I gives the observable quadrupole moment Uij 
including all terms up to the 2.5PN order. Therefore, by the previous reasoning, 
we can simply add to (5.10) of paper I the contribution of tails of tails at 3PN order 
and obtain the complete Uij to the 3PN order. We have 



U ii{ T R ) = Mg> + 2^ drM$\T R - r) 

c Jo 



ln Gro) 
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+ T2 
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+ OO 



dr 



M (3) .M^ 3) 

lyl a<i lyl j>a 



(Tr - r) - ^ a %M^ 



+ 2 



r+00 

/ drM^(T R -r) 
Jo 



( CT 



57 , f cr 



ln z — + — In — + 



\2r J 70 V 2r o 



+ ^ 



124627 
44100 

(4.13) 



The various terms are: at 1.5PN order, the dominant tail integral [32]; at 2.5PN 
order, the quadrupole-quadrupole terms [including in particular a non-local (memory) 
integral] and a quadrupole-dipole term (see paper I and references therein); and, 
at 3PN order, the tail of tail integral computed in this paper. The formula (4.13) 
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constitutes the main result of this paper, as it gives all the physical effects in the 
radiation field measured by a far away detector up to the 3PN order. 

Note that Uij, when expressed in terms of the intermediate moments Ml and 
Sl as in (4.13), shows a dependence on the (arbitrary) length scale r . Most of this 
dependence comes from our definition (4.6) of a radiative coordinate system, and thus 
can be removed by inserting Tr = t — r/c — (2GM j c 3 )ln(r /tq) back into (4.13), and 
expanding the result when c — > oo, keeping the necessary terms consistently. In doing 
so one finds that there remains a ro-dependent term at 3PN order, namely 



Uij = M\f - ^ In (j-^J (^P) Mf + terms independent of r . (4.14a) 

This term results simply from our use of the r -dependent formulas (3.4) and (3.11) in 
constructing the harmonic-coordinates metric. As we see from (4.14a), the dependence 
of Uij on ro (or rather r /c) is through the effective quadrupole moment 



This moment is exactly the one which appears in the near-zone expansion of the 
external metric, when taking into account the appearance of the dominant logarithm 
of c arising at the 3PN approximation. See the discussion in the Appendix of [29]. 
The appearance of this In c was pointed out by Anderson et al [44] . 



4.3. The energy flux 

We now investigate the total energy flux generated by the isolated source (its 
total gravitational luminosity) , and notably the non-local contributions therein. From 
(4.10) the luminosity in terms of the observable moments reads 



+oo 



c = V -9- I + „(!)„(!) , 41(1 + 2) (1) 

c 2t+i \(£-i)M(2£+l)U L L (£ — !)(•£+ 1)!(2£ + l)!!c 2 



1=2 



(4.15) 

The powers of 1 /c show notably that the moments Ul and Vl of higher multipolarities 
£ contribute to higher orders in the post-Newtonian expansion of C. 
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Let us gather the available information on non-local effects present in the U^s 
and Vl's. Rather, we consider the time-derivatives of the moments (U^ and V^), 
because these are the quantities of interest in (4.15). From (4.13) we write the time- 
derivative of the quadrupole Uij as 



UP = M? + 2™J*-*rM§> 



2r J 12 



H — g | instantaneous terms M a< iMj >a , and e a b<iMj >a Sb 

/GM\ 2 /-+ 00 , , 6) T 2 /cr\ 57, /cr\ 124627 
1 2 < — ) 1 ,' n U J + 70 U J + 44100" 







+ -^<j instantaneous terms M ab<l M j>ab , M ab M ijab , and e ab<i M J>ac S bc 



+ <>(£), (4-16) 

where we indicate only the index structure of the local (instantaneous) terms, but 
write m extenso all the non-local integrals. Note the important fact that the non- 
local (memory) integral present in at the 2.5PN order is a mere time anti- 
derivative [see (4.13)], and therefore becomes instantaneous when considering the 
time-derivative. We have added with respect to (4.13) the information that the 3.5PN 
term is instantaneous, exactly like the 2.5PN term. This follows from the dimensional 
argument used before. The 3.5PN term is such that 3(n — 1) + 2k = 7, therefore 
n = 2 and k = 2. Now we know [32] that the only non-local integrals at the quadratic 
order n = 2 are the tail integral which is purely of order 1.5PN, and the memory 
integral which contributes to the 2.5PN, 3.5PN and higher orders but is in the form of 
a simple time anti-derivative of an instantaneous functional of the moments Ml and 
Sl- This proves that the 3.5PN term in (4.16) is indeed instantaneous. Furthermore its 
multipole structure follows from i\ +£ 2 + s — 2 = 2k = 4. Note that the instantaneous 
terms in (4.16) (and other equations below) are instantaneous functionals not only 
of the moments Ml and Sl but also of the real source variables, i.e. when Ml and 
Sl are replaced by their explicit expressions as integrals over the source. Indeed the 
non-local integrals in Ml and Sl are not expected to arise before the 4PN order. 

Similarly we write the relevant higher-order multipole moments, but for them we 
need less accuracy than for the quadrupole. The results are 

25 



+00 



drM^ k 



In 



ct \ 97 
2n)J + 60 



+ ^-<J instantaneous terms M a<ij M k>a , e ab<i M jk>a S b 



and e ab<i M ja S k>b } + O 



(4.17a) 
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+ instantaneous terms M <j jMy > f + O ( — ) , 
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+ instantaneous terms M <i jS , / c> , and s ab< iMj a M k>b j + O — 

(4.17d) 

[The tails in the mass octupole and current quadrupole are from (5.8) in [35]. The mass 
2 4 -pole and current octupole are computed (including all the instantaneous terms) in 
(5.11) and (5.12) of paper I.] 

The separation made in the moments (4.16)-(4.17) between instantaneous and non- 
local terms yields a similar separation in the energy flux given by (4.15). Furthermore 
we introduce in the non-local part of £ a separation between the tail terms strictly 
speaking, a term involving the square of the tail, and the tail of tail term. Accordingly, 
we denote 



£ — £i n st + Aail + £(tail) 2 + Aail(tail) • (4-18) 

Quite evidently we include into the tail part of the flux, £taii> all the terms which are 
made of the cross products of the moments Ml and Sl and of the corresponding tail 
integrals at 1.5PN order. From (4.15)-(4.17) we obtain, neglecting terms at the 4PN 
order, 
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(4.19) 



The moments in front of each integrals depend on the current time T R . [For 
convenience we include in (4.19) the terms associated with the constants 11/12, 97/60, 
etc., though these terms are actually instantaneous. In fact these terms are given by 
some total time-derivatives and thus do not participate to the loss of energy in the 
source. In the case of binary systems moving on circular orbits, these terms are 
rigorously zero (see Section 5).] The expression (4.19) generalizes to 3.5PN order the 
expression (5.12) of [35]. 

Now the "(tail) 2 " contribution to the flux is given by the square of the tail integral 
at 1.5PN order, and therefore enters the energy flux at the same 3PN order as the 
contribution of tails of tails. In fact this contribution could be treated on the same 
footing as the "tail(tail)" contribution, but it will be clearer in Section 5 to investigate 
it separately. We have 
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(4.20) 



Physically (4.20) represents the energy flux due to the tail part of the wave, in 
situations where the tail can be separated from the other components of the field. 
In particular, after the passage of a burst of gravitational radiation (defined by the 
constancy of the quadrupole moment My before and after a certain interval of time), 
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the wave tail will be solely present in the radiation field, and therefore the total energy 
flux C will reduce in this case to £( ta ii) 2 - 

Next, the "tail(tail)" contribution to the flux involves the cross product of the 
quadrupole moment and of the tail of tail integral at 3PN order. It reads (neglecting 
4PN-order terms) 

\2r J 70 \2r J 44100 J \c 8 J J 

(4.21) 

Note that in contrast to (4.19) the instantaneous terms associated with the constants 
11/12 and 124627/44100 in (4.20) and (4.21) do contribute to the energy flux, even in 
the case of binary systems moving on circular orbits (see Section 5). 

Finally the instantaneous part of the flux, £i nst , is entirely defined, up to the 3.5PN 
order included, by the general formula (4.15) together with the previous definitions 
of Aaii, £(taii) 2 an d £taii(taii)- We do not write the full expression of £i ns t in terms of 
Ml and Sl because we do not consider it in the application to compact binaries in 
Section 5. It suffices for our purpose to recall that £i ns t not only is instantaneous in 
terms of the moments Ml and Sl, but also is instantaneous in terms of the real source 
parameters (neglecting 4PN-order terms), that is, in the case of compact binaries, of 
the orbital separation and relative velocity of the two bodies. 



5. Application to inspiraling compact binaries 

Last but not least we specialize the results to binary systems of compact objects 
(neutron stars or black holes). These systems, when the two objects spiral very rapidly 
toward each other in the phase just prior to the final coalescence (the orbital motion 
is highly relativistic in this phase), constitute the most interesting known source of 
gravitational waves to be observed by VIRGO and LIGO. The rate of inspiral is fixed 
by the total energy in the gravitational waves generated by the orbital motion, that is 
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by the binary's total gravitational luminosity £, which is therefore a crucial quantity 
to predict. We assume two non-spinning point masses (without internal structure) 
moving on an orbit which evolved for a sufficiently long time to have been circularized 
by the radiation reaction forces. For such (excellent) modelling of inspiraling compact 
binaries we compute Aaii, £(taii) 2 an d Aaii(taii) to the 3.5PN order included. 

To compute the tail part [given by (4.19)] we need the expressions of the multipole 
moments for circular compact binaries, to 2PN order for the mass quadrupole moment 
My, 1PN order for the mass octupole Mij k and current quadrupole Sij , and Newtonian 
order for Mj jfci and S^k (we need also the ADM mass M to 2PN order). These 
moments have been calculated in [45], and we simply report here their expressions: 
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M ljkl = umx <i:ikl> (l-3u)+o(^j , 
S ijk = urn e ab<i x jk>a v b (l -3u)+0 (J^j 



(5.16) 
(5.1c) 

(5.1d) 

(5.1e) 

(5-1/) 



The mass parameters are the total mass m = mi + 777-2, the mass difference 
5m = mi — m2, and the mass ratio v = 77717772/777 2 (satisfying < v < 1/4). The 
relative position and velocity of the two point-masses are denoted by x % = y\ — y\ and 
v l = dx l /dt. We use the post-Newtonian parameter 



Gm 
rc 



r = x 



(5.2) 



(r is the harmonic-coordinates distance between the two bodies). The time-derivatives 
of the multipole moments (5.1) are computed using the equations of motion, with 
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maximal 2PN-precision needed for the time-derivatives of the quadrupole moment. 
The equations of motion are 



dv % ( 1 

— = -ul PN x l + O^J , (5.3a) 
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1 + (-3 + i/) 7 + ( 6 + + v 2 ) 7 2 



(5.36) 



where the frequency W2PN is the orbital frequency of the exact circular periodic motion 
at the 2PN order (see [45]). Introducing instead of 7 the post-Newtonian parameter 
x = (GVmu^pn/c 3 ) 2 / 3 we have 

7 = *{l+(l-0*+(l-^y + O(* 3 )} • (5.3c) 

The computation of £taii proceeds like in Section VLB of [35]. Namely we insert into 
(4.19) the time-derivatives of the moments (5.1), evaluated both at the current time 
T R and at all anterior times T R — r. Then we work out all contractions of indices 
and obtain many integrals of logarithms of cr /2ri, where r\ is some constant such 
as T\ = r e -11 / 12 , multiplied by cosines of some multiples (n) of u^pnt- All these 
integrals are computed using the mathematical formula 



r+00 
Jo 



dylnye- Xy = ~(C + \nX) , (5.4) 
'0 ^ 

where A denotes the complex number A = 2ino;2PNri/c, and where C = 0.577.. is the 
Euler constant (see e.g. [47] p. 573). See Appendix A of [30] for the proof that this 
formula yields the correct result for inspiraling compact binaries modulo some error 
terms of order 0(c~ 5 lnc), falling in the present case into the un-controlled remainder 
of (4.19). In fact we need only the real part of (5.4) for the computation of Aaii, which 
leads [using ln(ia;) = Iniv + iir/2] to a term proportional to it and independent of r\. 
The result (extending the equation (6.18) in [33]) is 
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or, equivalently, in terms of the parameter x, 
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The tails contribute only to the half-integer post-Newtonian approximations 1.5PN, 
2.5PN, and 3.5PN. Now, we know that only the terms given by some non-local integrals 
can contribute to the half-integer post-Newtonian approximations. This follows from 
an argument presented in Sect. VI B of [33], which shows that the terms given by 
instantaneous functionals of the binary's relative position and velocity are zero for 
half-integer approximations in the energy flux for circular orbits (but only in this 
case). Thus we conclude that the terms computed in (5.5a) and (5.5b) represent 
the complete 1.5PN, 2.5PN, and 3.5PN approximations in C - no other contributions 
can come from £i ns t to these orders (and £( ta ii) 2 and £taii(taii) are purely of 3PN 
order). Being complete these approximations can thus be compared, in the test-mass 
limit v — > for one body, with the result of black-hole perturbation theory derived 
up to 4PN order by Tagoshi and Sasaki [50]. For the comparison we must use (5.5b) 
expressed in a coordinate-independent way by means of the parameter x. We find that 
there is perfect agreement, in the limit v — > 0, between (5.5b) and the corresponding 
terms in the equation (43) of [50] (see also the equation (3.1) of [51]). 

Turn now to £( ta ii) 2 defined by (4.20). The computation is essentially the same 
as for Aaii, but since we are considering a small 3PN effect the expression of the 
quadrupole moment at the Newtonian order is sufficient. Consistently we use the 
Newtonian equations of motion (with orbital frequency uj 2 = Gm/r s ). Using the 
mathematical formula (5.4), in which both the real and imaginary parts are now 
needed, we readily obtain 



-(tail) 2 
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C + ln 
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(5.6) 



Finally we compute £taii(taii) as given by (4.21). Again we need only the Newtonian 
quadrupole moment and Newtonian equations of motion, however a new ingredient is 
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necessary, which is a mathematical formula analogous to (5.4) but able to deal with 
the square of the logarithm. From [47] p. 574 the relevant formula is 



P + OO j 

J dy\n 2 ye- Xy = - 



n 
~6 



+ (C + lnA) ; 



(5.7) 



Using this formula for inspiraling compact binaries is justified in the same way as for 
the earlier formula (5.4). As a result we get 



--tail(tail) 
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(5.8) 



Both £(taii) 2 and Aaii(taii) have the same structure, namely that of a polynomial 
of the second degree in the combination C + ln(4u;ro/c). However, we see that the 
coefficients in front of the square of C + ln(4u;ro/c) in (5.6) and (5.8) are exactly 
opposite, and therefore that the sum of £( ta ii) 2 an d Aaii(taii) is actually merely linear 
in the combination C + ln(4a;ro/c). This fact is somewhat surprising because the 
terms involving the square of In a; are a priori allowed. Thus, to the 3PN order, the 
terms (lncu) 2 cancel out and there remains simply a term with In a; (and the associated 
lnc). The appearance of a In u at 3PN was first shown in this context by Tagoshi and 
Nakamura [39] . Recall that the general structure of the post-Newtonian expansion of 
the near-zone metric involves besides the regular powers of 1/c some arbitrary powers 
of lnc [16]. Similarly one expects that the general structure of C should involve when 
going to higher post-Newtonian approximations some arbitrary powers of In a;. [Note 
that up to the 5.5PN order C as computed in the limit v — > by Tanaka et al [50] is 
still linear in In a;.] 

Adding up (5.6) and (5.8) we obtain 
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(5.9) 



It is not yet possible to compare this result when v — > with the one obtained using 
the perturbation theory. Indeed (5.9) represents only a part of the complete 3PN term 
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in £, which should also take into account the (hard to get) 3PN contributions in £i nst , 
which are due to the 3PN relativistic corrections in the quadrupole moment (5.1a) 
and in the equations of motion (5.3). (The 3PN contributions in £i nst are in progress 
[56,57].) Nevertheless, we already recognize in (5.9) the same terms with it 2 and the 
combination C + ln(4a;) as in the result of Tagoshi and Sasaki [50]. [To the 3PN order 
one can replace in (5.9) the parameter 7 by x.] 

The 3PN term obtained in (5.9) depends on the arbitrary length scale tq. This 
is not a problem because the 3PN term in C- inst is expected also to depend on tq 
through in particular the explicit expression of the intermediate quadrupole moment 
My as a functional of the source's stress-energy tensor. The ro-dependent terms in 
both £(taii) 2 +taii(taii) an d Anst should cancel out, so that the physical energy flux C 
is indeed independent of r . We leave for future work [56,57] the check of the latter 
assertion. 



Appendix 1. Formulas to compute the cubic non-linearities 

In this Appendix we compute the finite part of the retarded integral of a source 
term with multipolarity £, radial dependence with k > 1, and containing a non-local 
integral whose kernel is a certain Legendre function Q m (x) [see (2.9)]. Thus, 



k ^ L = FP^qD^ 1 



/+00 
dxQ m (x)F(t - rx) 



(A.l) 



We tackle first the case k = 1. It can be checked in this case that the source 
behaves when r — > is such a way that the retarded integral (in its usual triple 
integral form) is convergent, so we can forget about the finite part procedure. For the 
present computation it is convenient to use the particular formula (D5) of Appendix 
D in [16], which gives 



t — r t+r-Z +00 

i, m ^L = ~ f <% [ 2 dw f dxQ m (x)F(£-(x-l)w) 
2r J.oo Ji-ipl Ji 



x P, 1 - 



(t - r - Q (t + r 
2rw 



i-2w) 



(A.2) 
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where Pi denotes the Legendre polynomial. Changing the variables (£,w) to the new 
variables (y, z) defined by £ — (x — l)w = t—ry and z = l — (t — r—^)(t+r—^ — 2w)/2rw, 
we obtain 



rriL 
'~2~.#i 



+00 



dx 



/+00 pi 
dyF(t-ry) J dz P e (z) 



x 



xy — z 



(A.3) 



A/(xy - z) 2 - (x 2 - l)(y 2 - 1) 

Next we integrate by parts the y-integral, introducing the anti-derivative F* -1 ) of F. 
After some manipulations we get 



f. r+00 p+00 
,m^L = -^J dyF^(t - ry) J dxQ m (x) 



x 



dx 



dz P t {z) 



1 y/{xy - z) 2 - (x 2 - l)(y 2 - 1) 



(A.4) 



The ^-integration can be performed explicitly thanks to the rather interesting 
mathematical formula 



dzP e (z) 



Pe(x)Qi(y) 
Pe(y)Qe{x) 



(l<x<y) 



(A.5) 



2 7_! ^( X y- Z )2-( X 2-l)(y2-l) (P e (y)Q t (x) (K y < x) . 

As this formula does not seem to appear in standard text-books of mathematical 
formulas such as [47] we present its proof at the end of this Appendix. 
With (A.5) we obtain 



i, m ^L = n L J^ dyF(-V(t-ry)^Qi(y) J\xQ m {x)^{x) 

+ P £ (y) J^°° dxQ m (x)^(x)^ . (A.6) 

The leading-order term at infinity (r — > 00, t — r = const) comes from the second term 
in the brackets, and we have 

i,m*L = — / drF^it-r-r) d xQ m (x)^(x) + o {r s ~ 2 ) (A.7) 

r JO Jl+r/r MX 
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(see the notation concerning the remainder in Section 3.2). In the case £ = m (the only 
one needed in this paper), the x-integral is obtained in closed- form, and we obtain 



i^l = ~J ^ drF^ 1 \t-r-r)[Q t (l+ T -)]\o{r £ - 2 ) . (A.8) 
Using the expansion (2.9c) of the Legendre function we further obtain 



tf L =- — / drF^it-r-r) 
° r Jo 



X 



2 I 



+ o(r £ - 2 ) . (A.9) 



This is the formula used in Section 3.2. 

Consider now the case k > 2, and restrict attention to the leading-order term at 
infinity. In this case one uses Lemma 7.2 in [16], whose hypothesis are satisfied with 
N = k — e > 2 — e. Thus the 1/r term in k,m^ 'l, which is deduced from (7.2) and 
(7.6) in [16], is 



(A.10) 



where u = t — r and 



G(u) = LTV-o / dsRl X) \^--l i8 



(A.lla) 
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/+oo 
dxQ m (x)F(s - (x - l)y) 



(A.116) 



By combining together these expressions, and introducing the variable z = l — (u—s)/y, 
we get 



1 f + °° f + °° ( y\ B 

G ( u ) = -7^pr Fp s=o J dxQ m (x) J dy I — I y 



-k+e+2 



i 

x / dz (z 2 - l) l F{u - (x - z)y) . 



(A.12) 
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Still this expression is not suitable for practical computations, and we must integrate 
by parts the y-integral in order to determine the finite part at B = 0. The all- 
integrated terms during the integrations by parts are zero by analytic continuation in 
B. We must distinguish two cases, 2 < k < £ + 2 and k > £ + 3. 

When 2 < k < £ + 2 we obtain after £ + 2 — k integrations by parts the local 
expression 



£\ 

where the coefficient k,m,otz reads 



GW(«)=- fcim a,F( fc - 3 )(«), (A.13) 



(£-k + 2)\ f + °° ^ M [\ (l-z 2 Y ,. _ 

k,ma e = 2e+l£l ^ dxQ m {x) J ^ dz j- _ ^_ fc+3 . (A.14) 

A more elegant form of this coefficient is obtained by using the representation of the 
Legendre function given by 

i r 1 (i-z 2 ) 1 

Q '^ = 2^L dz ¥^ (A ' 15) 

(which is equivalent to the other formulas (2.9), see e.g. [52]). Then k,m,cte reads 



/ f (z — x) 6 

dxQ m {x) J dz ^_ 3 j, Qi(z) ■ (A.16) 

The numerical values of the coefficient are computed from 



k 2 (-Y(k-2)l (2£ - 2k + 3 + 2i)\\ 



Q i\(k-i-2)\ (2£ + l + 2i)\\ 

/+oo 
dxQ m (x)Q e _ k+2+ 2i(x) , (A.17) 



where the remaining integrals are given by (see e.g. [47]) 



r+oo 
Jl 



, , (m-p)(m+p+l) L^J = 1 ) S=l i {m^p), 

dxQ m {x)Q p {x) = < r 2 L n J (A.18) 

1 ! ^-EUM (m = p). 

We thus obtain the coefficients reported in Tables 1 and 2 of Section 3. 

In the case k > £ + 3 we get after k — £ — 2 integrations by parts the non-local 
expression 
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where the coefficients k,mPe an d k,mle are given by 
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More elegant forms are 
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k,mPe = - J dxQ m (x) J dz ^_ 3 j, P i 
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1 f + °° f 1 (z-x) k ~ 3 

k,mlt = g J dxQ m (x) J dz—^——^-Pi(z) 



The cases k = 3 and £ = needed in this paper yield 
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dxQmix) 
r X to(to + 1) 
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(to > 1) 



(A.22a) 
(A.226) 



Finally the formula 7.132.3 in [47] where, however, a factor 2 2A M in the denominator 
should be corrected to a factor 2 2_/i , leads to 



3,m70 



18 T 12 



m +3m+3 



m(m+l)(m+2) \ m+1 

When m = 2 we find the values (3.13). 



2 v^m-l 1 
m — 1 Z^j = l j 



(m = 1) , 
(m > 2) . 



(A.22c) 



37 



Proof of the mathematical formula (A. 5 ). The author proved this formula by (i) 
verifying that the left-hand-side of the formula, namely 



_ 1 r 1 dzP e (z) 

= 2 y_ 1 _ z y _ p _ fyp _=fy ' 

is a particular solution of the Legendre equation both in variables x and y, and (ii) 
showing that this particular solution Ig(x, y) is necessarily equal to the solution given 
by the right-hand-side. To prove (ii) one must invoke the behaviour of the Legendre 
function at infinity, which is Qe(x) ~ l/x i+1 when x — > oo. 

Then a direct and more interesting proof of this formula was found by H. Sivak. 
This proof consists first of transforming Ii(x,y) by means of the formula 3.613.1 of 
Gradshteyn and Ryzhik [47] in order to obtain 

r / N 1 T r f 1 dzPAz) /A 

= / eft / i J — . (A.24) 

Z7T y ,/_! xy - ^ - a/ (x 2 - l){y 2 - 1) cost 

With the help of Neumann's formula (2.9a) for the Legendre function we get 

h(x, y) = - [ dt Q e (xy - ^ (x 2 - l)(y 2 - 1) cost) . (A.25) 
n Jo 

Using now the formula 8.795.2 of [47] we find (assuming for instance that x < y) 



h(x,y) = - / dt \P e (x)Q e (y) + 2y2(-) k P e - k (x)Q k e (y) coskt . (A.26) 
77 Jo L k=i ) 

This yields immediately the desired result 

I e (x,y) = P e (x)Q e (y) . (A.27) 
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Appendix 2. The harmonicity algorithm in the far zone 



The harmonicity algorithm is a mean to add to the finite part of the retarded 
integral of the source term, 



u 



af3 = FPb^D- 1 



ro 



B 



A a/3 



(B.1) 



an homogeneous solution v af3 of the wave equation which is such that the sum u a ^ '+v al3 
is divergence- free (and thus satisfies the harmonic-gauge condition). A particular 
algorithm is defined by (4.12)-(4.13) in [16]; a slightly different one by (2.11)-(2.12) in 
paper I. Here we follow paper I. 

The harmonicity algorithm to order 1/r when r — > oo, t — r = const is as follows. 
The 1/r term of the divergence of u al3 ', 



W = V(n,t-r)+o(l 
is decomposed into multipole moments according to 



(B.2) 



U° = Y,n L W L , (B.3a) 
U l = ^2n iL X L + ^2{n L - 1 Y iL - 1 + e iab n aL - 1 Z bL - 1 } . (B.36) 



£>0 £>1 



Comparing this decomposition with the definition (2.11) in paper I, we obtain some 
relations between the tensors Wl, Xl, etc., and the time-derivatives of the tensors 
A L , B L: etc., in this definition. Next we obtain the 1/r term of w a/3 , 

v aP = -V Q/3 (n,t-r) + (J-^J , (B.4) 

by applying the formula (2.12) in paper I. Re-expressing the decomposition in terms 
of the tensors Wl, Xl, etc., we obtain some formulas for the time-derivative of V al3 , 



dV 00 

— = -W-n a [W a + Y a - 3X a ] , (B.5a) 
— — = -Yi + 3X, - e iah n a Z h + ^ n L -iW iL -i , (B.56) 

£>2 
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dt 

e>2 



5ijX + -25 ij n L - 1 X L _ 1 + 6n L _ 2 (iX j)L _ 2 

+ n L - 2 [WijL-2 - 3X ijL - 2 + YijL-2] 

+ 2n aL - 2 £ab(iZj)bL-2^ • (B.5c) 



These are the formulas needed in Section 4.1. 



Acknowledgment s 

The author would like to thank Horacio Sivak for his proof of the mathematical formula 
(A. 5) and his permission of presenting it at the end of Appendix A. 



40 



References 



[1] Blanchet L Quadrupole-quadrupole gravitational waves (referred to as paper I) 

[2] Foures-Bruhat F 1952 Acta Math. 88 141 

[3] Bruhat Y 1964 Ann. Math. Pura Appl. 64 191 

[4] Kundt W and Newman E T 1968 J. Math. Phys. 9 2193 

[5] McLenaghan R G 1969 Proc. Camb. Philos. Soc. 65 139 

[6] Friedlander F G 1975 The wave equation on a curved space-time (Cambridge: 

Cambridge University Press) 
[7] Waylen P C 1978 Proc. R. Soc. London A 362 233 

[8] Choquet-Bruhat F, Christodoulou D and Francaviglia M 1979 Ann. Inst. Henri 
Poincare A 31 399 

[9] Carminati J and McLenaghan R G 1986 Ann. Inst. Henri Poincare (Physique 
Theorique) A 44 115 
[10] Bonnor W B 1959 Philos. Trans. R. Soc. London A 251 233 
[11] Bertotti B and Plebanski J 1960 Ann. Phys. (N.Y.) 11 169 
[12] Thorne K S and Kovacs S J 1975 Astrophys. J. 200 245 
[13] Crowley R J and Thorne K S 1977 Astrophys. J. 215 624 
[14] Thorne K S 1980 Rev. Mod. Phys. 52 299 

[15] Bel L, Damour T, Deruelle N, Ibanez J and Martin J 1981 Gen. Relativ. Gravit. 
13 963 

[16] Blanchet L and Damour T 1986 Philos. Trans. R. Soc. London A 320 379 

[17] Blanchet L 1987 Proc. R. Soc. Lond. A 409 383 

[18] Peters P C 1966 Phys. Rev. 146 938 

[19] Price R H 1972 Phys. Rev. D 5 2419 

[20] Price R H 1972 Phys. Rev. D 5 2439 

[21] Bardeen J M and Press W H 1973 J. Math. Phys. 14 7 

[22] Schmidt B G and Stewart J M 1979 Proc. R. Soc. London A 367 503 

[23] Porrill J and Stewart J M 1981 Proc. R. Soc. London A 376 451 

[24] Leaver E W 1986 Phys. Rev. 34 384 

[25] Bonnor W B and Rotenberg M A 1966 Proc. R. Soc. London A 289 247 
[26] Couch W E, Torrence R J, Janis A I and Newman E T 1968 J. Math. Phys. 9 
484 



41 



[27] Hunter A J and Rotenberg M A 1969 J. Phys. A 2 34 

[28] Bonnor W B 1974 in Ondes et radiations gravitationnelles (CNRS, Paris, 1974) 
p. 73 

[29] Blanchet L and Damour T 1988 Phys. Rev. D 37 1410 

[30] Blanchet L and Schafer G 1993 Class. Quantum Grav. 10 2699 

[31] Blanchet L and Sathyaprakash B S 1995 Phys. Rev. Letters 74 1067 

[32] Blanchet L and Damour T 1992 Phys. Rev. D 46 4304 

[33] Blanchet L 1995 Phys. Rev. D 51 2559 

[34] Will C M and Wiseman A G 1996 Phys. Rev. D 54 4813 

[35] Blanchet L 1996 Phys. Rev. D 54 1417 

[36] Cutler C, Apostolatos T A, Bildsten L, Finn L S, Flanagan E E, Kennefick D, 
Markovic D M, Ori A, Poisson E, Sussman G J and Thorne K S 1993 Phys. 
Rev. Lett. 70 2984 

[37] Cutler C, Finn L S, Poisson E and Sussmann G J 1993 Phys. Rev. D 47 1511 
[38] Tanaka T, Shibata M, Sasaki M, Tagoshi H and Nakamura T 1993 Prog. Theor. 
Phys. 90 65 

[39] Tagoshi H and Nakamura T 1994 Phys. Rev. D 49 4016 
[40] Poisson E 1995 Phys. Rev. D 52 5719 

[41] Our notation is the following: signature h ++; greek indices =0,1,2,3; 

latin indices =1,2,3; g = det (g^); rj^ = j]^ u = flat metric = diag (- 
1,1,1,1); r = |x| = (x\ + x\ + x^) 1 / 2 ; n l = Hi = x 1 /r; <9; = d/dx l ; 
n L = ul = ni 1 rii 2 ...ni e and = d^d^ . . .di e , where L = i^.-.ii is 
a multi-index with i indices; n^-i = . . . n a L-i = n a %_i, etc. . . ; 

fiL and £>l are the (symmetric) and trace-free (STF) parts of ul and <9l, 
also denoted by n<L>, <9<l>; the superscript (n) denotes n time derivatives; 
T( a p) = \{T a i3 + Tp a ) and T (ij) = \{T i5 + T Si ) 

[42] Fock V A 1959 Theory of Space, Time and Gravitation Pergamon, London 

[43] Madore J 1970 Ann. Inst. H. Poincare 12 285 and 365 

[44] Anderson J L, Kates R E, Kegeles L S and Madonna R G 1982 Phys. Rev. D 
25 2038 

[45] Blanchet L, Damour T and Iyer B R 1995 Phys. Rev. D 51 5360 
[47] Gradshteyn I S and Ryzhik I M 1980 Tabie of Integrals, Series and Products 
(Academic Press, Inc.) 

42 



[48] Poisson E 1993 Phys. Rev. D 47 1497 

[49] Sasaki M 1994 Prog. Theor. Phys. 92 17 

[50] Tagoshi H and Sasaki M 1994 Prog. Theor. Phys. 92 745 

[51] Tanaka T, Tagoshi H and Sasaki M 1996 Prog. Theor. Phys. 96 1087 

[52] Whittaker E T and Watson G N 1990 A course of Modern Analysis reprinted 

edition (Cambridge U. Press) 

[53] Leonard S W and Poisson E Phys. Rev. D in press |gr-qc/9705014 

[54] Blanchet L and Damour T 1989 Ann. Inst. H. Poincare 50 377 

[55] Damour T and Iyer B R 1991 Ann. Inst. H. Poincare 54 155 

[56] Blanchet L, Iyer B R and Joguet B in preparation 

[57] Blanchet L, Faye G, Iyer B R, Joguet B and Ponsot B in progress 

[58] Damour T, Iyer B R and Sathyaprakash B S gr-qc/9708034 



43 



Table 1: Values of the coefficients k,m,oti for k — 2. 



Table 2: Values of the coefficients k, m oze for k — 3. 
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